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1 Introduction 



One of the most interesting branches of field theory is devoted to spatially 
extended topological solitons in 3+1 dimensional space-time such as domain 
walls or vortices. By now it is clear that these solitons, mathematically 
represented by certain nontrivial solutions of classical field equations, play 
an important role in many phenomena in condensed matter physics P], EL H], 
particle physics and perhaps also in cosmology |J. Dynamics of the 
extended solitons is extremely rich, but also it poses rather hard problems 
because of its strongly nonlinear character. In contradistinction to the case 
of solitons in one- dimensional space, there are no methods of constructing 
generic exact domain wall or vortex solutions. Till now most theoretical 
results have been obtained with the help of numerical methods, see, e.g. 
I@j IU> eventhough there are certain remarkable exact analytical results like 
the discovery of travelling waves on a straightlinear vortex ||, or of the 90° 
scattering of straightlinear parallel vortices JT0 . 

Analytical approaches to evolution of a generic curved vortex/domain 
wall are scarce. In most of them one attempts to reduce the dynamics to 
much simpler dynamics of a string/membrane. This line of thought has been 
initiated in the seminal paper jnj, continued in JT2"] , F3 and later on in [14- 



21]. It has been found already in [11], |T^] that in the case of local vortices, e.g. 
in the Abelian Higgs model, one may approximately describe motion of the 
vortex using the simple classical Nambu-Goto string model as the leading 
approximation. Dynamics of internal degrees of freedom of the vortex is 
completely neglected in this approximation. In a reference frame co-moving 
with the vortex all fields have by assumption the same form as for a static, 
straightlinear vortex. In other words, the field structure of the vortex in the 
co-moving reference frame is kept frozen. The case of domain walls is strictly 
analogous — here one has the Nambu-Goto membrane instead of the string. 

In the case of a global vortex a long range field is present. In the leading 
approximation the collective dynamics of the vortex is described in terms of 



the Nambu-Goto string coupled to a massless 2-form gauge field [p_^ 



Our present paper is devoted to the problem how to compute corrections 
to the collective dynamics which are due to the internal dynamics of the vor- 
tex/domain wall. Specifically, we further develop the approach proposed in 
[p3|| . In the literature one can also find other attempts to calculate evolution 
of the curved vortex/domain wall beyond the leading approximation. Let us 
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review them briefly. 

In the most popular approach [14-21] the string/membrane is identified 
with the locus of zeros of a scalar field, that is with the core. Influence 
of the internal dynamics on the collective motion is taken into account by 
corrections to the Nambu-Goto action for the core. The total action defines 
the so called effective classical string /membrane model. It turns out that 
such effective models are rather unpleasant. Equations of motion contain 
higher derivatives with respect to time and consequently they have unphysical 
runaway type solutions. Moreover, it seems that the effective models should 
be nonlocal, see [H] where caveats concerning this approach are presented 
and Section 4 of the present paper. 

In another approach [2S| one derives directly from field equations approx- 
imate equations of motion for certain averages, e.g. for an average radius in 
the case of a ring-like vortex, in an analogy with the Ehrenfest's approach to 
evolution of a quantum mechanical wave packet in an external potential. It 
is clear that the knowledge of time evolution of such averages is not sufficient 
in order to recover the time evolution of the fields of the vortex/domain wall. 

In the polynomial approximation the string/membrane is identified 
with the core, but it gives only a part of the collective dynamics. The other 
collective degrees of freedom are represented by coefficients of certain polyno- 
mials which approximate the fields inside the vortex/domain wall. They can 
be regarded as scalar fields defined on the worldsheet of the string/membrane. 
Here one avoids shortcomings of the previous two approaches. On the other 
hand, the accuracy of this approximation has been checked only in the sim- 
plest cases of static straightlinear vortices and planar domain walls. More- 
over, it is likely that this approach should be generalised to include radiation 
of massive fields from the extended solitons. Such radiation could dump os- 
cillations of the width which in the present formulation of the method can 
persist indefinitely (this problem is currently under investigation). 

The approach proposed in |23| is based on the expansion in the width 
for the vortex/domain wall solution of the field equations. It is constructed 
along the lines of Hilbert-Chapman-Enskog method for singularly perturbed 
evolution equations | 2~7fl . In this approach certain consistency conditions 



play a crucial role. Co- moving string/membrane is of the Nambu-Goto type 
in all orders, but it does not coincide with the core. In addition to the 
string/membrane there are also scalar fields defined on the worldsheet of 
the co-moving string/membrane. We avoid the unpleasant effective classical 
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models with higher derivatives. In contradistinction to the case of polynomial 
approximation we have here physically meaningful expansion parameter, that 
is the width defined as the inverse mass (M _1 ) of the scalar field. More 
precisely, the expansion parameters are given by dimensionless ratios of the 
width to local curvature radia of the string/membrane worldsheet. In the 
case of a domain wall the terms in the expansion have been calculated up 
to the order M~ 2 . These results have been compared with purely numerical 
solutions for cylindrical and spherical domain walls |28| — it has turned out 



that the approximate analytical solution is astonishingly accurate. 

Expansion in the width was invoked also in the earlier papers, e.g., in 
|L4], [L5|, [17], |2l| where a version of it was used on intermediate stages of 
derivations of the effective string/membrane models. However, in those pa- 
pers the consistency conditions were not taken into account. This is a crucial 
difference because in our version of the expansion just these conditions make 
the calculation of the effective classical string/ membrane action superfluous. 
More detailed comparison of our version of the expansion with the other 
ones can be found in the first of papers [p]|1 . See also Section 4 of the present 
paper. 

The construction of the expansion given in [23] has two main elements: a 
special co-moving coordinate system jl2] introduced to maintain the Lorentz 
invariance order by order in the expansion, and the consistency conditions. 
While the use of the co-moving coordinates could be regarded as a rather 
elegant and relatively simple way to secure the Lorentz invariance, the con- 
sistency conditions appeared in [23] as formal, purely mathematical obser- 
vations. In spite of the crucial role they played in determining the vor- 
tex/domain wall solution, their field theoretical meaning was hidden. More- 
over, the calculations were rather cumbersome and even for the domain wall, 
which is simpler than the vortex, we were able to obtain only the first and 
the second order corrections. We found that the core did not coincide with 
the co-moving string/membrane, and that the co-moving string/membrane 
obeyed Nambu-Goto equation in all orders, but we did not derive any equa- 
tion of motion for the core. Such equation would be particularly useful 
for comparison with the effective string/membrane approach in which the 
string/membrane is identified with the core. 

In the present paper we concentrate on the case of a domain wall. Analo- 
gous results can be obtained also for a vortex but the necessary calculations 
are much longer. We reformulate the derivation of the expansion in such a 
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way that the (2+l)-dimensional scalar fields defined on the worldsheet of the 



co-moving membrane, which have been introduced in |23| for purely mathe- 
matical reasons, acquire a physical interpretation — they are identified with 
perturbative contributions to the values of the original scalar field on the 
membrane. We also find a field theoretical interpretation of the consistency 
conditions: they coincide with certain Euler-Lagrange equation. Moreover, 
derivations of the higher order terms in the expansion are significantly sim- 
pler in the new formulation. Now we can consider corrections in the next 
two orders, that is in the M~ 3 and M -4 orders. Order by order in the 1/M 
expansion we express profile of the domain wall by the collective degrees 
of freedom given by the co- moving membrane and the (2+l)-dimensional 
scalar fields. Finally, we derive an approximate equation of motion for the 
core. It turns out to be different from equations obtained in the effective 
action approaches. In particular, it does not contain higher derivatives. 

The main conclusion of our paper is that the collective dynamics of the do- 
main wall beyond the leading approximation involves the (2+l)-dimensional 
scalar fields defined on the worldsheet of the Nambu-Goto membrane. Elim- 
ination of those fields would lead to a nonlocal membrane model for the 
collective dynamics. 

The plan of our paper is as follows. In Section 2 we present Euler- 
Lagrange equations and we discuss certain important consequences of them. 
In Section 3 we consider the extrinsic curvature corrections to the domain 
wall solution up to the M -4 order. The equation of motion for the core is 
derived in Section 4. Section 5 is devoted to a discussion of our results. In 
Appendix A we have collected formulas relevant for the transformation to 
the co-moving coordinates. In Appendix B we present solutions of certain 
linear ordinary differential equation. In Appendix C we list functions and 
constants which appear in Section 3. 

2 The consistency conditions as Euler — La- 
grange equation 

We consider a domain wall in the model defined by the following Lagrangian 
£ = -- % „d^d»$ - ^ - ^) 2 , (1) 
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where $ is a single real scalar field, (r/ AtI/ )=diag(-l, 1,1,1) is the space-time 
metric, and A, M are positive constants. Two vacuum values of $ are equal 
to ±M/2 V / A. The mass of the corresponding scalar particle is equal to M. 

The domain wall appears when a solution $ of the Euler-Lagrange equa- 
tion obtained from the Lagrangian (|l|) smoothly interpolates between the two 
vacuum values. Then, at each instant of time the field $ vanishes somewhere 
in the interior of the domain wall. The locus of these zeros is assumed to be a 
smooth connected surface S in the space Q. It is called the core of the domain 
wall. The transverse width of such a domain wall is of the order M~ l and 
energy density is exponentially localised around the core. The world-volume 
£ of the core is a 3-dimensional manifold embedded in Minkowski space-time. 
The core plays an essential role in the effective membrane models mentioned 
in the Introduction. 

In our approach we use another smooth connected surface S attached to 
the domain wall. In general it differs from the core, except at the initial 
instant of time when it coincides with the core by assumption. We shall 
obtain Nambu-Goto equation for S, so this surface can be regarded as the 
Nambu-Goto type relativistic membrane co-moving with the domain wall. 
Let us stress that this co-moving membrane is merely an auxiliary mathe- 
matical notion. It is introduced in order to define the co-moving coordinate 
system. The world- volume of S is denoted by S. We shall parametrise it as 
follows 

E 3 (Y»)(u a ) = (r,Y i (u a )). (2) 

We use here the notation (w a ) a =o,i,2 = (t, a 1 , a 2 ), where r coincides with the 
laboratory frame time x , while a 1 , a 2 parametrise the co-moving membrane 
S at each instant of time. The index i — 1, 2, 3 refers to the spatial compo- 
nents of the four-vector. The points of the co-moving membrane S at the 
instant r are given by (F*)(r, a 1 , a 2 ). The coordinate system (r, a 1 , a 2 , £) 
co-moving with the domain wall is defined by the formula 

i^yyj + ^V), (3) 

where x M are Cartesian laboratory frame coordinates in Minkowski space- 
time, and (n M ) is a normalised space-like four-vector orthogonal to £ in the 

1 We use slightly different notation than in [^3) — the present one is more convenient, 
especially in Section 4. 
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covariant sense, 



0, n M n M 



where Y% = dY^/du a . The three four- vectors Y a are tangent to E. The 
definition (3) implies that £ and u a are Lorentz scalars. In the co-moving 
coordinates the co-moving membrane is described by the simple Lorentz in- 
variant condition £ = 0. For points lying on S the parameter r coincides 
with the laboratory time x°, but for £ ^ in general r is not equal to x°. 
In the co-moving coordinates Lagrangian (1) has the following form 



C 



M 4 r 4 

[«9 S 0«90 S + (0 2 - l) 2 + j^G ab d a ^ 



32A 



(4) 



where d s = d/ds,d a = d/du a , and G ab is given in Appendix A. The field $ 
and the coordinate £ have been expressed by dimensionless and s, 



According to papers p3f the core is shifted with respect to the co-moving 



membrane, that is the scalar field has non-zero values on the membrane. We 
take this fact as the starting point for the new derivation of the expansion 
in the width. The idea consists in extracting from the scalar field its compo- 
nent living on the co-moving membrane and treating it separately from the 
remaining part of the scalar field. To achieve this we use the identity 

( p(s,u a ) = B(u a )^s) + X (s,u a ), (5) 

where 

B(u a ) I (j)(0,u a ) (6) 
is the component of the scalar field living on the co-moving membrane, and 

X d ± ^u a ) - B{u a )4,{s) (7) 

is the remaining part. The auxiliary, fixed function ip{s) depends on the 
variable s only. It is smooth, concentrated around s — 0, and 

^(0) = 1. (8) 
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It follows that 

x (0,u a ) = 0. (9) 

We shall see that the best choice for ip(s) is given by formula (22) below. 
Formulas (5-7) can be regarded as an invertible change of variables 

<f>(s,u a )^(B(u a ), X (s,u a )) 

in the configuration space of the scalar field. Therefore, it is legitimate to 
use formula (5) in Lagrangian (4) and to derive Euler-Lagrange equations 
by taking independent variations of B(u a ) and x- The variation 5x has to 
respect the condition (9), hence 

S X (0,u a )=0. 

Because of this condition, variation of the action functional 




with respect to x gives Euler-Lagrange equation in the regions s < and 
s > 0. It has the following form 

-L d a W=gh&*db(Bi> + x)\ (io) 

M z \J—g 

+ l -d s [hd s {B^ + X )\ + KB^ + x )[l - (B4> + xf] = 0, 

where h, g and G ab are given in Appendix A. At s = there is no Euler- 
Lagrange equation corresponding to the variation 5x- Instead, we have the 
condition (9). Equation (10) should be solved in the both regions separately, 
with (9) regarded as a part of boundary conditions for x- To complete the 
boundary conditions we also specify the behaviour of x f° r I CI niuch larger 
than the characteristic length 1/M, that is for \s\ ^> 1. In our model (1) 
the expected behaviour of the domain wall field $ for large |£| is given by 
an exponential approach to the vacuum values. Therefore, we shall seek a 
solution such that x is exponentially close to +1 for s 3> 1 , while for s <^ — 1 
it is exponentially close to —1. 

At this stage of our considerations Eq.(10) should not be extrapolated to 
s = 0. For example, the l.h.s. of it could have a <5(s)-type singularity. It 



7 



would occur if x was smooth for s > and for s < but had a spike at 
s = 0. 

In addition to Eq.(10) we also have the Euler-Lagrange equation corre- 
sponding to variations of B{u a ). This equation has the following form 



— I 

M 2 J 



ds -^d a \^-ghG ah d b {B^ + x)] ■*!> (11) 
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-\fds hd s ^d s {B^ + x) 

+ Jds W{Bii + x) [l-W + xf 



= 0. 



Here and in the following we use /rfsasa shorthand for the definite integral 

Comparing equations (10) and (11) one might think that they are not 
independent because it seems that multiplying Eq.(10) by ip and integrating 
the result over s with the help of integration by parts we obtain Eq.(ll). 
This argument is false, namely it ignores the above mentioned possibility 
that the l.h.s. of Eq.(10) might have the 5-type singularity at s = 0. On the 
other hand, if the singularity at s = is absent then indeed Eq.(ll) does 
follow from Eq.(10). 

Actually, just because Eqs.(10) and (11) are independent we can prove 
that the solution x does not have the spike at s = 0. To this end, we 
multiply (10) by ip and integrate over s in the intervals (oo, e], [— e, — oo) 
with a positive e which approaches zero. In this manner we avoid the not- 
excluded-yet singularity at s = 0. The l.h.s. of the resulting formula is 
compared with the l.h.s. of Eq.(ll): in the limit e — > they differ by the 
term ^ 

--*phd s x\s=o+ + -?phd s x\s=o-, 

which has to vanish because the r.h.s.'s of Eqs.(10), (11) vanish. In conclu- 
sion, 

lim dsx = lim d s x, 

s— >0+ s^O— 

that is d s x is continuous at s = 0. Hence, Eq.(10) and Eq.(ll) together imply 
that the spike at s = is absent. It follows that Eq.(10) is obeyed by x & l so 
at s = 0. To summarize, Eq.(ll) and the statement that the singularity at 
s = is absent are equivalent. 
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Let us now solve Eq.(lO) in the leading approximation which is obtained 
by putting 1/M = 0. The equation is then reduced to 

^ S V°) + 0( O )[1-(0(°)) 2 ]=O, (12) 

where 

0(°) = jB ( o V + x (o) . 

In the 1/M expansion this is the only nonlinear equation we have to solve. 
Mathematically, it coincides with a well-known equation for a planar domain 
wall. It has the following particular solutions 

<f) ( ±\s,u a ) = ±tanh(s - b{u a )), (13) 

where b(u a ) is an arbitrary function of the indicated variables. In the fol- 
lowing we shall take the + sign (the — sign corresponds to an anti-domain 
wall). The functions b(u a ) can be removed from the solution (13) by a suit- 
able choice of the co-moving membrane. The new membrane is obtained by 
shifting the points of the worldsheet E of the original membrane along the 
direction normal to S by the distance b(u a ). Because the normal direction is 
given locally by the n{u a ) four-vector, this shift is equivalent to the change 
(£,w a ) — > (£',w a ), £' = ^ — b(u a ) of the co-moving coordinates. Thus, we may 
take as the zeroth order domain wall solution in the co-moving frame 

(o) =tanhs, (14) 

without any loss of generality. This solution together with conditions (8), 
(9) gives 

B {0) = 0, x (0) =tanhs. (15) 

The solution (14) in the co- moving coordinates does not determine the 
field in the laboratory frame because we do not know yet the position of the 
co-moving membrane with respect to the laboratory frame. Equations (10), 
(11) should also yield an equation for the co- moving membrane, otherwise 
they would not form the complete set of evolution equations for the domain 
wall. In fact, we shall see that the first order terms in Eq.(10) imply the 
Nambu-Goto equation for the membrane. 

The expansion in the width in the present formulation has the form 

X (s, u a ) = tanhs + ^X (1) (s, u a ) + -^X i2) (s, u a ) + t^X (3) (s, u a ) + (16) 
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B(u a ) = + ^5 {2) K) + ^ (3 V) + (17) 

where we have taken into account the zeroth order results (15). The expan- 
sion parameter is 1/M and not 1/M 2 because 1/M in the first power appears 
in the h and G ab functions after passing to the s variable, see formulas in 
Appendix A. In order to obey the condition (9) and to ensure the proper 
asymptotics of x at large \s\ we assume that for n > 1 

X (n) (0,u a ) = 0, lim x (n) = 0. (18) 

Inserting the perturbative Ansatz (16), (17) in Eqs.(10) and (11), expanding 
the l.h.s.'s of them in powers of 1/M, and equating to zero coefficients in 
front of the powers of 1/M we obtain a sequence of linear, inhomogeneous 
equations for x^(s, u a ), B^ n \u a ) with n > 1 . 

The first order terms in Eq.(10) give the following equation 

L X {1) = -B^lAj} + K a a d s x {0 \ (19) 

where 

L + 1 - 3( X <°>) 2 , 

and x^ is given by the second of formulas (15). General solution of the 
equation of the form (19) is presented in Appendix B. 

The most important point in our derivation of the expansion in the width 
is the observation that the operator L has a zero-mode, that is the normal- 
izable solution 

^ o(s) = c^k 

of the homogeneous equation 

Lip = 0. (20) 

Notice that = d s x^ — this means that the zero-mode i[)q is related to 
the translational invariance of Eq.(12) under s — > s + const. The presence 
of the zero-mode implies the consistency condition: we multiply Eq.(19) by 
■00 and integrate over s. With the help of integration by parts we show that 
/ ipoLx^ and / ipoLip vanish because of (20). Finally, we obtain the following 
condition 

K a a jds Ms)dsX { °\s)=0 
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which is equivalent to 

K a a = 0. (21) 

This condition coincides with the well-known Nambu-Goto equation. It de- 
termines the motion of the co-moving membrane, that is the functions Y l (u a ), 
i=l,2,3, once initial data are fixed. When we know these functions we can 
calculate the extrinsic curvature coefficients and the metric g a ^. Review 
of properties of relativistic Nambu-Goto membranes can be found in, e.g., 

Due to Nambu-Goto equation (21) one term on the r.h.s. of Eq.(19) 
vanishes. The resulting equation 

L X (1) = -B {1) L^ 

has the following solution obeying the boundary conditions (18) 

It can be obtained from the formulas given in Appendix B, but there is 
a shorter way: notice that + B^'ip obeys the homogeneous equation 
(20) which has the general solution C {u a )ip$(s) + D(u a )ipi(s); the boundary 
conditions (18) imply that D = 0, C = Now, we recall that the function 
ip(s) is an auxiliary mathematical object, like the co-moving membrane in our 
approach, and we may adjust it in order to simplify formulas. The obvious 
choice is 

= Ms) = -4l-> ( 22 ) 
cosh s 

because then the first order correction has the simplest possible form 



X 



il) {s,u a )=0. (23) 



In the remaining part of our paper we will use ip given by formula (22). 
Notice that vanishing x does not mean that the first order correction to 
the original field <fi also vanishes — there is the first order contribution equal 
to B^i/jq/M. It does not vanish on the co-moving membrane that is at s = 0. 

Equations (10), (11) in the first order do not give any restriction on the 
function B' 1 ', In the next Section we shall see that an equation for B^ 
(Eq.(32) below) follows from the third order terms in Eq.(ll). This situation 
is typical for singular perturbation theories of which the 1/M expansion is 
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an example - - higher order equations imply restrictions (the consistency 
conditions) for the lower order contributions [^?J . 

With the choice (22) for if), Eq.(ll) expanded in the powers of 1/M gives 
equations for B^ n > coinciding with the consistency conditions found in [p3[| . 
This follows from the fact that both Eq.(ll) and the consistency conditions 
are obtained by multiplying expanded with respect to 1/M Eq.(lO) by the 
zero-mode ip$ and next integrating over s. Euler-Lagrange equation (11) can 
be regarded as the generating equation for the consistency conditions. 



3 The higher order corrections 

At this point we have the complete set of equations determining the evolution 
of the domain wall in the 1/M expansion. Each of Eqs.(lO), (11), (21) 
describes a different aspect of the dynamics of the curved domain wall. We 
shall see that expanded in the powers of 1/M Eq.(lO) determines dependence 
of x on s, that is on the distance from the worldsheet E of the co-moving 
membrane along the perpendicular direction n(u a ) at each point Y(u a ) of 
the worldsheet. Because the term Bip in formula (5) has explicit dependence 
on s, see formula (22), we may say that Eq.(lO) for \ fixes the transverse 
profile of the domain wall. 

Equation (11) determines the B^{u a } functions, which can be regarded 
as a (2+l)-dimensional scalar fields defined on the worldsheet S and having 
nontrivial nonlinear dynamics, see Eqs.(32), (35) below. The extrinsic cur- 
vature of £ acts as an external source for these fields. The fields B'"' 
can propagate along E. One may regard this effect as causal propagation of 
deformations which are introduced by the extrinsic curvature. 

Finally, Nambu-Goto equation (21) for the co-moving membrane deter- 
mines the evolution of the shape of the domain wall. 

Equations (21) and (11) are of the evolution type in the 1/M expansion 
- we have to specify initial data for them, otherwise their solutions are not 
unique. Equations for the perturbative contributions x are °f different 
type — in order to ensure uniqueness of their solutions it is sufficient to 
adopt the boundary conditions (18). The initial data for B{u a ) and Y % {u a ) 
follow from initial data for the original field <fi. From such data for <f> we know 



the initial position and velocity of the core. As in papers [23| we assume that 



at the initial instant tq the co-moving membrane and the core have the same 
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position and velocity. Hence, 



initial data for the membrane = initial data for the core. 

Using formula (5) one can show that then 

B^(r ,a\a 2 ) = 0, d T B^(r Q , a\ a 2 ) = 0. (24) 

Our perturbative scheme with the initial conditions (24) gives the basic curved 
domain wall. In Section 5 we discuss more general solutions obtained by 
adopting more general than (24) initial data for the B^> fields. 

In order to find the domain wall solution one should first solve the col- 
lective dynamics, that is to compute evolution of the co-moving membrane 
and of the B field. The profile \ °f the domain wall is found in the next step 
from formulas (16) above and (27) below. In our perturbative scheme the 
profile of the domain wall can not be chosen arbitrarily even at the initial 
time — it is fixed uniquely once the initial data for the membrane and for the 
B field are given. Evolution of the core can be determined afterwards, from 
the explicit expression for the scalar field <fi. This will be done in Section 4. 

Now we shall show that order by order in the 1/M expansion Eq.(10) de- 
termines the dependence of the field x on s. x wm De explicitly expressed by 
the perturbative contributions to the B field and by the geometric character- 
istics g a b, of the worldsheet of the co-moving membrane. Therefore, the 
dynamics of the curved domain wall is reduced to the Nambu-Goto dynamics 
of the co-moving membrane and to nonlinear dynamics of the B^ fields. A 
convenient starting point for this calculation is obtained by rewriting Eq.(10) 
in a different form. We substitute 



where 



with 



and 



X = X (0) +X, h = l + h, G ab = g ab + G ab , 



h = —h {2) + — /i (3) 
= -2s 2 K\Kl = - 8 -s*KZK b c K:, 

Qab = J_ G (l)ab + _L G (2)ab + Qf M -3\ 

M M 2 v ' 
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with 



Q(l)ab = AsK ab q(2) 



ah 



■Yls 2 K a K ca . 



We also use Eqs.(12), (20) and formulas (15). In the end Eq.(10) is written 
in the following form 



M 2 



--d a 



g (l + h){g ab + G ab )d b {B^ + x) -hL X (25) 



■dshtyo + dsX + Bd s ij ) + (1 + h)(3 X {0) + Bij + x)W + x? ■ 



Expanded in the powers of 1/M it gives equations of the type considered in 
Appendix B, that is 

L X {n) =f {n) . (26) 

The source term is determined by the lower order terms in x an d B, 
namely contains x^ with k < n — 2 and with I < n — 1. Solution 
of Eq.(26) is given by the formula 



X (n) (s) = JdxG(s,x)f n \x), 



(27) 



where the Green's function G(s, x) can be found in Appendix B. This solution 
obeys the boundary conditions (18). 
Let us list /W with n=2,3,4: 



/( 2 ) = 2# ^ + 3 X (0 Vo 2 (5 (1) ) 2 , 



(28) 



/( 3 ) = As 2 ^ K^K b K c a - 2^ D (3) 5 (1) 

+fi 3) mZK b a BW + /f (a )( S (i))3 + 6x ^ Q B^B^, 



r(4) 



where 



; (4) = _ 2D (3) (x (2) +^(2)) _4 S ^ 1 d^y/^K^B^) 



'-0 



(29) 



(30) 



-l^ww - ±0.Zi< 2 >(0.iM*< 2 > + 9 sX (2) - /> (2 Vo) 
+3^ 2 (5( 1 )) 2 (x (2) + ^ (2) ) + 3 X <°W(5< 2 >) 2 + (X (2) ) 2 + 2^oX (3) 5 (1) ], 
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and 

/ 2 (4) = 6^X (0) 5 (1) 5 (3) . 

Here 

□ (3) ^ -?=d a (V^g ab d b ) 



-9 

is the three-dimensional d'Alembertian on the world-volume S of the co- 
moving membrane. The functions /^(s) can be found in Appendix C. The 
formulas for J*- 3,4 ** will become a little bit simpler because we shall show that 
vanishes. 

From formula (27) we obtain the explicit dependence of x^ 011 s - The 
integrals over x which appear on the r.h.s. of that formula after we substi- 
tute expressions (28-30) are rather simple. Most of them can be evaluated 
analytically and the rest numerically by a computer algebra system. From 
the remark at the end of Appendix B and from formula (28) we see that 
is an odd function of s. The correction x^ is & n even and x^ an °dd 
function of s provided that vanishes. 

Equations of motion for the (2+l)-dimensional fields B^{u a ) are ob- 
tained by expanding Eq.(ll). It is convenient to rewrite that equation in the 
following form 



— I 

M 2 J 



ds *Po-^d a k/=£(l + h)(g ab + G ab )d b (B^ + x)] (31) 



9 

1 f 

+- J ds d s h(i>l + xd s ip + B^ d s ^o) 

- j ds ^(1 + h)(3x (0) + Bij + x)W + x) 2 = 0. 

It has been obtained with the help of integration by parts. We have also used 
Eq.(20). From Eq.(31) we see that the first nontrivial equation appears in the 
order 1/M 3 . It has the form of non-linear, inhomogeneous (2+l)-dimensional 
wave equation for B^(u a ) regarded as a field defined on £ 

^□(3)5(1) + CxK a h K h a B^ - c 2 {B^ f = d K^K b K c a . (32) 

The constants q, do are listed in Appendix C. 

The equation for B^ follows from the 1/M 4 terms in Eq.(31). It has the 
form 

□ (3)5(2) + \ K a K b aB {2) _ C A( B mf B i2) = Q (33) 

2 c 
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This equation is homogeneous with respect to B^ 2 \ Because also the initial 
data (24) are homogeneous, B^ vanishes 

B^(u a ) = 0. (34) 

The 1/M 5 terms give the equation for B@\ It can be written in the form 

coDWb® + Cl K a b K b a B^ - 3c 2 (B^) 2 B^ = , (35) 

where the source term does not contain B@\ It is given by the following 
formula 

j (5) = - Ids -^d a iV=g(g a %x {3) ^ + G (1)a W 2 tyo (36) 

J v 9 
+h^g ab d b B^m + G^ ab d b B^m)} 

-\ J ds d.Md.h {S) X {2) + d s h^x {3) ) + \jds ^)(B«) S 

+ 3 - J ds M^xW + ^x (3) (£ (1) ) 2 + Mx (2) ) 2 b ( v 

+i,l x ^h^\B^f + 2^ (2) X (0) X (2) 5 (1) ] 
+3jds^GUP]B^, 



where 



G[A\s) = JdxG( 8 ,z)fl A \x). 



The function B^ is not present in formula (36) because we have taken into 
account the result (34). 

Finally, from Eq.(31) in the 1/M 6 order we obtain a homogeneous wave 
equation for the function B^ and therefore, similarly as in the case of the 
function B^ with Eq.(33), we conclude that 

B (4) K) = 0. (37) 
Let us summarize our results. We have found that 
= X (0) (s) + iB«(u%(«) + -^X {2) (s,u a ) (38) 
^\u a )Ms) + -^X {3) (s,u a ) + 4lX (4) (*X) + 0(M~ 5 ), 



M 3 v /ruv ' M 3 v ' M 4 ' 
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where x^°\^Po are given by formulas (15), (22), are given by formula (27), 
and are solutions of Eqs.(32), (35) with the initial data (24). Before 

solving Eqs.(32), (35) one should first solve Nambu-Goto equation (21) for 
the co-moving membrane. This is necessary in order to determine g ab and 
K ab which appear in these equations, and also in order to find the explicit 
form of the transformation (3) which relates the co-moving coordinates to 
the Cartesian laboratory coordinates. 



4 Equation of motion for the core 

In the description of the collective dynamics of the domain wall we could use 
the core instead of the co-moving membrane. This might even seem a better 
choice because the definition of the core as the locus of zeros of the scalar 
field is very simple and it directly refers to the domain wall solution. In 
this Section we derive the equation of motion for the core in our approach. 
We shall see that this equation is more complicated than the Nambu-Goto 
equation for the co-moving membrane and that it involves the B^ fields. 

The worldsheet £ of the core S is related to the worldsheet £ of the 
co-moving membrane by the formula 

£ 3 Y»{u a ) = Y»{u a ) + (39) 

where s is determined from the equation 

(f)(s,u a ) = 0. (40) 

Taking for <f> the expansion (38) we find that 

5(1) 

s = -— + 0(M- 3 ). (41) 

We parametrise the core by the same coordinates (u a ) as for the co-moving 
membrane. The induced metric tensor and the extrinsic curvature for £ are 
defined by the formulas 

(jab = Y a Y b , K ab = -n, a Y b = nY ab , (42) 
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where the normal to E 

nK) = n(«°) - + 0(M~ 4 ) 

has been determined from the condition 

nY a = for a = 0,1, 2. 

We have also used the Nambu-Goto equation K° = 0. Simple calculations 
give 

~9ab = 9ab + -wB^K ah + 0(M" 4 ), (43) 
K ab = K ab + ^BWK ac K c b - -^V b B^ + 0(M~ 4 ), (44) 

where 

is the standard covariant derivative of B$ with respect to the metric g a b on 
the co-moving membrane. The relations inverse to (43), (44) have the form 

9ab = ~9ab ~ j^B^K ah + 0(M- A ), (45) 
K ab = K ab - ^B^K ac K c b + -^V b B^ + 0(M" 4 ), (46) 

where the tilda over V means that now we use the metric g ab . Using formulas 
(45), (46) we find that the equation K% = for the co-moving membrane 
implies the following equation for the core 

K = ~^- 2 B^m b a - ^ Q(3) £ (1) + 0{M^). (47) 
Using Eq.(32) we eliminate n (3) B {1) . Then 

K = -^(3 - C j)B^ktkl - ^^(5 (1) ) 3 (48) 

2 do -KK b c k: + o(M~*). 



M 2 c ' 

This is the final form of the equation of motion for the core. 

The terms on the r.h.s. of Eq.(48) can be regarded as corrections to the 
Nambu-Goto equation K% = 0. They vanish if KlK h c K c a = because then 
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also B^ = as it follows from Eq.(32) with the initial conditions (24). In 
general, Eq.(48) has to be considered together with Eq.(32) for B^ in which 
we may replace K ab by K ab . It is clear that the preferred by us description of 
the collective dynamics in terms of the co-moving membrane of the Nambu- 
Goto type instead of the core is much simpler. 

In the literature one can find many attempts to represent the collective 
dynamics of the domain wall by the core only. In our approach this would 
amount to expressing the B^ field in Eq.(48) by the extrinsic curvatures. 
In principle this is possible because Eq.(32) with the initial data (24) gives 
a one-to-one relation between B^ and K b K h c K c a . However, the presence on 
the l.h.s. of Eq.(32) of the operator acting on B^ has the consequence 
that in general B^ depends on K b K h c K c a in a nonlocal manner. Hence, the 
resulting self-contained equation for the core will be nonlocal too. 

5 Discussion 

1. The (2+l)-dimensional fields B^ n \u a ) are an important component in our 
version of the expansion in the width. The leading term B^ 1 ' /M does not 
vanish if KlK h c K c a ^ 0. Using formulas from the second of papers |26| one can 
show that KlK h c K c a is proportional to the product of the mean and Gaussian 
curvatures in a local rest frame of an infinitesimal piece of the co-moving 
membrane. It is also shown in that paper that the evolution governed by the 
Nambu-Goto equation is sensitive only to the mean curvature in the local 
rest frame. Therefore, the co-moving membrane, which obeys the Nambu- 
Goto equation to all orders, can not account for effects which are due to the 
Gaussian curvature. This is the reason why for the complete local description 
of the collective dynamics of the domain wall we need the B field too. 

2. There is an important assumption we have tacitly made: that the 
derivatives Xa-i-^a are °f the or der 1/M n . It is not satisfied, for example, 
if x an d B contain modes oscillating with a frequency ~ M which give posi- 
tive powers of M upon differentiation with respect to u a . If such oscillating 
components were present the counting of powers of 1/M would no longer be 
so straightforward as in Sections 2 and 3. The assumption excludes radia- 
tion modes as well as massive excitations of the domain wall. Therefore, the 
approximate solution we obtain gives what we may call the basic curved do- 
main wall. To obtain more general domain wall solutions one would have to 
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change appropriately the approximation scheme. Actually, the fact that such 
particular radiationless unexcited curved domain wall exists is a prediction 
coming from the 1/M expansion. The expansion yields domain walls of con- 
crete transverse profile — the dependence on s is explicit in the approximate 
solution we construct even at the initial instant of time. Once we choose the 
initial position and velocity of points of the membrane the dependence of 
the scalar field on the variable s at the initial time is given by formulas (16), 
(27). This unique profile is characteristic for the basic curved domain wall. 

3. The approximate domain wall solution we have obtained can be gen- 
eralised by relaxing the initial conditions (24) for the B^ fields. These 
homogeneous initial data forbid appearence of freely propagating along the 
co-moving membrane waves of the flW fields — only the fields gener- 
ated by the extrinsic curvature as the source can appear. Suppose however 
that we admit such free waves, that is that the initial data for B^ are cho- 
sen arbitrarily 0. Now B^ and B^ do not have to vanish. Equations (32) 
and (33) remain unchanged, but in Eq.(35) new terms containing B^ will 
appear. Let us also choose certain initial data for Nambu-Goto equation 
(21) for the co-moving membrane. Then, substitution of the corresponding 
unique solutions of Eqs.(21), (32), (33) and of (the modified) Eq.(35) in for- 
mula (38) gives an approximate solution (f) of the domain wall type. These 
more general solutions are not excluded by the assumption formulated in the 
preceding paragraph because the characteristic frequency of the B^ fields is 
not related to the mass parameter M. From Eqs.(32), (33), (35) we see that 
the B^ fields have the effective mass 2 ~ K^K b a . 

Actually, a simple particular case of solutions analogous to the generalised 
ones we are now discussing has been reported in paper JHJ]. In that paper 
certain small amplitude excitations of a straightlinear vortex in the Abelian 
Higgs model have been considered. Such excitations are represented by fields 



defined on the axis of the straightlinear vortex. They have been called in [30 



2 Of course, we have to respect the usual limitations imposed by consistency of the 
perturbative approach: 

r(«) r(«+i) . , r(") 

M n M M™ 

These bounds ensure that order of magnitude of the perturbative contributions, in principle 
given by a power of K a b/M, is not changed because of too large values of the fields 
at the initial time or later. 
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the zero-mode fields. They obey a linear, massless wave equation. It turns 
out that our field gives the domain wall counterpart of the zero-mode 



field of Ref.Bn] if we take a planar domain wall for which K ab = 0, and if 
we assume that the field is so small that one may neglect the nonlinear 
term on the l.h.s. of Eq.(32). 
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7 Appendix A. The co-moving coordinates 

The extrinsic curvature coefficients K ab and induced metrics g ab on E are 
defined by the following formulas: 

J^-ab — 'V ,a&> Mob — 1 )a J H,bi 

where a, b = 0, 1, 2. The covariant metric tensor in the new coordinates has 
the following form 

G a b 

1 



[Gap] : 

where a, /3 — 0, 1, 2, 3; a = 3 corresponds to the £ coordinate; and 



G ab = N ac g cd N db , N ac I g ac - ^K ac . 
Thus, = 1, G^ a = 0. Straightforward computation gives 



where as usual g — det[g a b], G = det[G a p], and 

u a ) = 1 - £Kl + \e{K a a Kl - K h a Kt) - \eKlK h c Kl. 

Also g can depend on u a . For raising and lowering the Latin indices of the 
extrinsic curvature coefficients we use the induced metric tensors g ab , g ao . 
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The inverse metric tensor G a @ is given by the formula 

[G al3 \ ■ 



G" 





1 



where 



G ab = ( N -i^ac gcd ( N -ijdb_ 
Explicit formula for (iV -1 ) ac has the following form: 



(AT 



-l\ac 



l {g ac [l - iK b b + \e{K b b Kj - K*K b d )} 
+ i{l-iK h b )K ac + eK a d K dc ). 



This is just the matrix inverse to [N ab \. It has the upper indices by definition. 

In general, the coordinates (u a ,^) are defined locally, in a vicinity of 
the world-volume £ of the co-moving membrane. Roughly speaking, the 
allowed range of the £ coordinate is determined by the smaller of the two 
main curvature radia of the membrane in a local rest frame. We assume that 
this curvature radius is sufficiently large so that on the outside of the region 
of validity of the co-moving coordinates there are only exponential tails of 
the domain wall, that is the field is exponentially close to one of the two 
vacuum solutions. More detailed discussion of the region of validity of the 
co-moving coordinates can be found in the first of papers 



8 Appendix B. Solutions of equation L\ = / 

Let us consider the equation 

^ 2 x+(l-3tanh 2 S )x = /(s). 

The presence of the zero mode implies that it has solutions only if f(s) obeys 
the condition 

Jds Ms)f(s)=0, 

where / ds is our shorthand for J_™. General solution of our equation can 
be obtained with the help of the standard procedure consisting of finding 
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two independent solutions of the corresponding homogeneous equation and 



constructing a suitable Green's function, see, e.g., [31|. The two linearly- 



independent solutions of the homogeneous equation have the form 

^o(s) 



2 „ ' 



cosh s 

■01 (s) = - sinh(2s) + - tanh s + - 2 



8 cosh s 
As the Green's function we take 

G(s, x) = 2 [V>i(s)V>o(aO - Ms)M x )] ( s ~ x ) + 2^ o (s)^i(x)0(-x), 

where B is the step function. This Green's function vanishes at s — 0. The 
general solution of our equation has the form 



X = aipo + j3ipi + I dxG(s, x)f(x), 

where a, (3 do not depend on s. If we require that x obeys the conditions 
(18) we have to put a — (3 — 0. In this case, using the given above formula 
for the Green's function we can write the solution x i n the following form 

/s ps 
dx 0o (x) f(x) - 2^o (s) / dxip x (x) f (x) . 
-oo JO 

This specific formula for x( s ) ^ s used in Section 3. One can easily check that 
if / is even (odd) function then x is even (odd) function too. 



9 Appendix C. The list of functions and con- 
stants 

It is convenient to introduce the notation 

G[f](s) = J dxG(s,x)f(x). 

The functions appearing in formula (29) for J® have the form 
A (3) ( S ) = 2sd s i> + 12x (0 VoG[x^ ], 
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/ 2 (3) ( S ) = ^ + 18 X (0 VoG[x {0 Vo 2 ]- 

The constants present in Eqs.(32), (33), (35) are given by the following for- 
mulas: 

d = 2 J ds s 2 ?pQ, c = J ds ipo, 
ci = ~ J ds ^oA (3) , 

° 2 = \ I ds ^ ^ (3) ' 
c 3 = lJds ^ + 18 J ds ^X (0) G[x (0 Vo 2 ]- 
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